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The Australia Telescope 20 GHz (AT20G) Survey: analysis 
of the extragalactic source sample 
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■ ABSTRACT 

The Australia Telescope 20 GHz (AT20G) survey is a blind survey of the whole South- 
ern sky at 20 GHz with follow-up observations at 4.8, 8.6, and 20 GHz carried out with 
the Australia Telescope Compact Array (ATCA) from 2004 to 2008. In this paper we 
present an analysis of radio spectral properties in total intensity and polarisation, sizes, 
optical identifications, and redshifts of the sample of the 5808 extragalactic sources 
in the survey catalogue of confirmed sources over 6.1 sr in the Southern sky (i.e. the 
k> , whole Southern sky excluding the strip at Galactic latitude |6| < 1.5°). 

5_j ■ The sample has a flux density limit of 40 mJy. Completeness has been measured 

as a function of scan region and flux density. Averaging over the whole survey area 
the follow-up survey is 78% complete above 50 mJy and 93% complete above 100 
mJy. 3332 sources with declination S < —15° have good quality almost simultaneous 
observations at 4.8, 8.6, and 20 GHz. The spectral analysis shows that the sample 
is dominated by flat-spectrum sources, with 69% having spectral index a|° 6 > —0.5 
(S oc v a ). The fraction of flat-spectrum sources decreases from 81% for S20 GHz > 500 
mJy, to 60% for £20 GHz < 100 mJy. There is also a clear spectral steepening at 
higher frequencies with the median a decreasing from —0.16 between 4.8 and 8.6 GHz 
to -0.28 between 8.6 and 20 GHz. 

Simultaneous observations in polarisation are available for all the sources at all 
the frequencies. 768 sources have a good quality detection of polarised flux density 
at 20 GHz; 467 of them were also detected in polarisation at 4.8 and/or at 8.6 GHz 
so that it has been possible to compare the spectral behaviour in total intensity and 
polarisation. We have found that the polarised fraction increases slightly with fre- 
quency and decreases with flux density. The spectral indices in total intensity and in 
polarisation are, on average, close to each other, but we also found several sources for 
which the spectral shape of the polarised emission is substantially different from the 
spectral shape in total intensity. The correlation between the spectral indices in total 
intensity and in polarisation is weaker for flat-spectrum sources. 

Cross matches and comparisons have been made with other catalogues at lower 
radio frequencies, and in the optical, X-ray and 7-ray bands. Redshift estimates are 
available for 825 sources. 
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1 INTRODUCTION 

The properties of radio source populations at high 
( ^ 10 GHz) radio frequency are poorly known due to the 
fact that large-area high-frequency surveys are very time- 
consuming with high sensitivity ground-based diffraction 
limited telescopes. The characterization of this population 
of radio sources is of great interest since it is dominated by 
the young flat-spectrum AGNs which are expected to have a 
strong influence on the evolution of the early Universe. Fur- 
thermore, it is crucial for the interpretation of data collected 
by experiments imaging the Cosmic Microwave Background 
(CMB). Fluctuations due to radio sources are the main con- 
taminants of the CMB signal on scales smaller than 30 ar- 
cmin (De Zotti et al. 1999, Toffolatti et al. 2005); they need 
therefore to be carefully subtracted to avoid biasing the esti- 
mates of cosmological parameters. Because of the complex- 
ity and variety of radio spectra, this subtraction cannot be 
properly done relying on extrapolations from low frequency 
surveys (Massardi et al. 2009). 

This has motivated large observational efforts in the 
last several years to conduct surveys at high radio frequency 
(e.g., Waldram et al. 2003, 2009; Muchovej et al. 2009; see 
De Zotti et al. 2010 for a review). The Australia Telescope 
Compact Array (ATCA)Q 20 GHz (AT20G) survey, carried 
out from 2004 to 2008, covered all the Southern sky (6.1 sr) 
to a flux density limit of ~ 40 mjy (Murphy et al. 2010). 

A Pilot Survey (Ricci et al. 2004a; Sadler et al. 2006) 
at 18.5 GHz was carried out in 2002 and 2003. It detected 
173 sources in the declination range —60° to —70° down to 
100 mjy and was used to optimize the observational tech- 
niques for the full AT20G survey. 

The analysis of the 320 sources of the AT20G Bright 
Source Sample (S20GHZ > 0.50 Jy) at declination below —15° 
was presented by Massardi et al. (2008). In this paper we 
extend the analysis to the full sample of confirmed extra- 
galactic sources from the Murphy et al. (2010) catalogue. The 
outline of the paper is the following. § [2] contains the de- 
scription of the sample. In § we analyze the complete- 
ness and derive source counts for the whole sample as well 
as for the steep- and flat-spectrum populations separately. 
The spectral properties in total intensity and polarisation 
are discussed in § [4] and [5] respectively. Optical identifica- 
tions, redshifts, and cross correlations with surveys in other 
wavebands are briefly presented and discussed in § [5] More 
detailed analyses are deferred to subsequent papers. The re- 
sults are summarized in § (7) 



2 SAMPLE DEFINITION AND PROPERTIES 

2.1 Observations and sample selection 

The AT20G survey is a blind survey using a wide band (8 
GHz) analogue correlator connected (operating in the fre- 
quency range 16-24 GHz) to three of the ATCA 22m an- 
tennas. The telescope scanned at 15° per minute in declina- 
tion along the meridian to cover 10° to 15° wide declination 
strips, with 10 mjy rms noise. The whole Southern sky was 
completely covered in four observing epochs between 2004 
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and 2007. A comprehensive description of the instrumen- 
tal setup, the scan strategy, the map-making and source ex- 
traction techniques will be presented in a forthcoming paper 
(Hancock et al., in preparation). 

As described in Murphy et al. (2010), candidate sources 
from the blind survey scans were re-observed in a compact 
hybrid configuration with the ATCA digital correlator with 
two 128 MHz bands centered at 18752 MHz and 21056 MHz 
and two orthogonal linear polarisations for each band, dur- 
ing several observing epochs between 2004 and 2008. These 
follow-up observations were used to confirm real detections 
and measure accurate flux densities and positions. The com- 
bination of the two close bands are considered as a single 
256 MHz wide band centered at 19904 MHz, which is the 
reference frequency for our '20 GHz' observations. The main 
follow-up was performed soon after the blind scans, and few 
further runs were done to recover bad weather observations 
or improve the sample follow-up completeness. Almost si- 
multaneously (i.e. within a couple of weeks) from any 20 
GHz follow-up, detected sources were re-observed with an 
East-West configuration with two 128 MHz bands centered 
at 4.8 and 8.64 GHz (hereafter referred as '5' and '8' GHz) 
for spectral analysis. A data-analysis pipeline was built to 
produce the final catalogues. Above S = —15° the beams are 
very elongated for East-West ATCA arrays used for the 5 
and 8 GHz observations so we did not observe the sources in 
this region at frequencies below 20 GHz. The primary beam 
FWHM is 2.4, 5.5, and 9.9 arcmin at 20, 8, and 5 GHz, 
respectively and the angular resolution for the follow-up is 
typically 10 arcsec at 20 and 5 GHz and 5 arcsec at 8 GHz. 
The details of follow-up observations and the data reduc- 
tion procedures, together with the catalogue of the surveyed 
sources were presented in Murphy et al. (2010). 

The AT20G survey sample of followed-up sources con- 
sists of 5890 sources over 6.1 sr of the Southern sky at 20 
GHz (i.e. the whole Southern hemisphere with a Galactic cut 
for |6| < 1.5°). A map of the distribution of these sources 
is in Figure 8 of Murphy et al. (2010). Optical counterparts 
were obtained from the NASA Extragalactic Database (see 
§[6]) and in the SIMBAD database: they have been used to 
identify the Galactic objects in the sample. The exclusion 
of the \b\ < 1.5° region in fact does not guarantee that all 
sources in the Murphy et al. (2010) catalogue are extragalac- 
tic. The 65 sources identified as PNe were found by looking 
for flat spectrum sources that were slightly extended, se- 
lected according to what discussed in i]2.3l A full analysis 
including this data will be published separately. The 6 HII 
regions were identified through a search of the literature. 
Finally, the catalogue includes 82 sources (1.4%) flagged 
as Galactic or Magellanic Cloud HII regions or planetary 
nebulae, which we removed from our sample, leaving 5808 
sources. We kept the other sources in the Magellanic Cloud 
region. The sample may still contain some Galactic sources, 
but we believe that the residual contamination should be 
small because the relatively high resolution 20 GHz beam 
preferentially selects compact sources (< 30 arcsec) and 
Galactic sources with \b\ > 1.5°) are rarely that small. 

The final sample of sources used for the following anal- 
ysis is composed of 5808 extragalactic sources. 
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2.2 Flux density measurements 

Flux densities of point-like sources were calculated using the 
triple product technique according to which the flux density 
of a point source is the geometric average of the visibil- 
ity amplitudes in a baseline closure triangle. The flux den- 
sity uncertainty is a combination of errors due to calibration 
and to the noise level (for further details see Murphy et al. 
2010). This technique is very robust to the effects of phase 
de-correlation. We measured both the triple correlation and 
the peak amplitude in the image. These agreed very well in 
good weather conditions. In periods of poor phase stability 
the scatter between multiple observations using image am- 
plitudes increased while the scatter in the triple correlations 
was unaffected. Hence, the triple product technique was used 
to get accurate flux density measurements for almost all the 
compact sources. Only 327 sources have the 'poor quality' 
flux density flag set, due to low quality data. 

The triple product technique, however, is not well- 
suited for extended or multiple-component sources. We have 
identified the extended sources by combining automatic 
methods and visual inspections of both images and visi- 
bility functions as described by Murphy et al. (2010). The 
automatic extended source identification method includes 
the comparison of flux density measured on short and long 
baselines, and the analysis of any significant deviation of the 
phase closure from zero which is expected for a point source. 
Sources with significant flux density from regions larger than 
5-10 arcsec are flagged as extended using these methods. 
Our sample contains 308 sources in this category; thus they 
constitute 5% of the total sample. For 246 of these sources 
the integrated flux density was obtained from the visibility 
amplitude measured on the shortest spacing. However, this 
method can fail for sources with multiple components sepa- 
rated by more than about 45 arcsec, for which cases we don't 
have close enough short spacings to get a reliable integrated 
flux density. For these remaining sources the triple product 
is used instead as a rough estimate of the flux density of the 
strongest component and these sources are flagged as 'poor 
quality' flux densities. Multiple widely separated compact 
components (e.g. hot spots in lobes) may be listed as dis- 
tinct sources and may be misidentified (see §[6j) . Sources with 
structure within the 2.4 arcmin 20 GHz primary beam are 
flagged as extended. The position in the catalogue for these 
sources corresponds to the peak in the image which is usu- 
ally the flat-spectrum core in this high frequency catalogue, 
and for these sources the optical identification is correct. 
However, for a very small number of extended sources the 
strongest peak could be a hotspot in the lobes or jet and 
the identification might be unclear. These extended sources 
and refined identifications will be discussed in a future paper 
(Mahony et al. in preparation). 

We must stress that the survey and the source detec- 
tion technique is optimized for point sources, and there is 
bias against extended sources with angular sizes larger than 
about 45 arcsec. The 30m spacing used for the blind survey 
underestimates the amplitude by 50% and the survey is in- 
creasingly incomplete for more extended sources. Integrated 
flux densities at 20GHz have been obtained for nine very 
extended extragalactic sources (Burke-Spolaor et al. 2009). 
The flux densities determined for the seven sources included 
in AT20G survey were added to the final catalogue. No in- 



tegrated flux densities at other frequencies are available for 
these sources at the same epoch; so they are not included in 
the spectral analysis. 

Sources found to be extended at 20 GHz have been as- 
sumed to also be extended at the lower frequencies and the 
flux density was estimated from the shortest baselines. 

2.3 6 km visibility data 

Data were collected with the '6 km' antenna which is used 
in conjunction with the hybrid arrays used for the ATCA 
during the follow-up observations. This provided 5 baselines 
from 4.3 to 4.5 km for 5542 (94%) of the sources in the 
sample. The visibility of an interferometric measurement of 
a source is the Fourier transform of the brightness distri- 
bution of the source and decreases with increasing angular 
size, from 1.0 corresponding to a point source. The visibility 
is determined by the ratio of the average scalar amplitude 
of the 5 long (~ 4.5 km) baselines to the average scalar am- 
plitude of the 10 short (30 - 214 m) baselines. Given the 
typical errors in these 4.5 km visibilities, any source with 
measured visibility < 0.85 is considered to be significantly 
resolved. This corresponds to a size of 0.15 arcsec for a Gaus- 
sian brightness distribution or 0.23 arcsec for double peaks. 
We refer to these sources as 'not-compact'. This information 
was used to separate the compact and not-compact sources 
discussed in the following sections. A full discussion of these 
high resolution observations will be given in Chhetri et al. 
(in preparation). 

Note that the limiting size for compact sources defined 
above (~ 0.2 arcsec) corresponds to a much smaller angular 
size than in the discussion of extended sources in the pre- 
vious section which is based on the observations with the 
compact hybrid array alone. The definition of 'not-compact' 
based on the 6 km antenna visibility corresponds to sizes 
of more than 0.2 arcsec and contains 19% of the 20 GHz 
population. As shown in §[4l this size limit divides the com- 
pact and extended source populations very cleanly and cor- 
responds closely to the traditional classification into flat- 
and steep-spectrum sources. 

2.4 Polarisation 

As described in Murphy et al. (2010), the polarised flux 
density has been measured simultaneously with the total 
intensity measurements at all frequencies, and polarisation 
values at any frequency are given if: 

• the polarised flux density is at least three time larger 
than its rms error; 

• it is larger than 6 mjy; 

• the fractional polarisation is larger than 1 per cent. 

All the other cases have been considered non detections and 
we used as an upper limit to the polarised flux density the 
maximum between: a) 3 times the rms error on the polarised 
flux density; b) 6 mjy; and c) 1 per cent of the total flux 
density. 

No estimation of integrated polarised flux densities was 
attempted for extended sources and these are not included in 
the polarisation analysis. This results in some bias, however 
to include them we would need to use the peak polarisa- 
tion as given in the Murphy et al. (2010) catalogue which 
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Figure 1. Orthographic projection of the Southern hemisphere. Colors represent the completeness above 50, 70 and 100 mjy respectively, 
from left to right. RA=0 is towards the top of the figure and increase in the anticlockwise direction. The grey region corresponds to 
the Galactic plane region (|6| < 1.5°) and other regions not considered in the completeness analysis. The region with 8 > —15° and 
16h < RA < 18h has been discarded because the follow-up survey was seriously affected by bad weather. Also the reliability estimation 
techniques fail close to the South Pole and for this reason we discarded the <5 < —80° region in the completeness analysis and source 
counts. 



would result in an even worse bias because it will always 
be higher than the integrated polarisation which is what is 
measured for the unresolved sources. The only exceptions 
are the integrated polarisations for the sample of bright ex- 
tended sources included in Burke-Spolaor et al. (2009). 



3 COMPLETENESS, RELIABILITY, AND 
SOURCE COUNTS 

3.1 Blind scan and follow-up completeness 

The follow-up observations were carried out at the position 
of candidate objects detected on the survey blind scans. 
Candidate sources were prioritized for follow-up observa- 
tions in order of decreasing flux density to optimize sample 
completeness at high flux density levels. However, this fact 
constrained the flux density limit to that which could be 
reached in the allocated follow-up time given the variable 
weather conditions. As a result, different regions of the sky, 
observed during different follow-up campaigns have slightly 
different flux density completeness limits. The final sample 
completeness is determined by both the effectiveness of the 
survey source detection techniques (i.e. by the completeness 
of the candidate source list) and by the quality (i.e. observ- 
ing strategy, weather conditions and sensitivity level) of the 
follow-up observations. 

As discussed in Hancock et al. (2010, in preparation), 
and Hancock (2010) the completeness of the blind scan sam- 
ple was assessed by inserting simulated sources into the data 
stream at the raw (uncalibrated) data level, and following 
these sources through all the calibration and source detec- 
tion steps. With this procedure Hancock (2010) found that 
the blind scan sample is 97% complete above 50 mjy and 
93% complete above 40 mjy. However, this candidate source 
list is only 50% reliable in some regions of the sky at its 
lower flux density limits. To obtain a catalogue with essen- 
tially 100% reliability we needed the follow-up observations 



but as discussed above, the follow-up observing strategy and 
data quality was the major limitation to the final sample 
completeness. To quantify the completeness as a function of 
flux density and region of the sky we can compare the list of 
confirmed sources in the follow-up catalogue with the blind 
scan survey source list whose completeness is known. We 
could not use this technique directly on the follow-up cat- 
alogue since it is affected by the weather conditions at the 
moment of the follow-up observations and by the efficiency 
of the detection algorithm that improved over time. How- 
ever, we could apply this method to the blind scan survey 
because all the data was re-reduced with a single consistent 
algorithm over the entire survey. Obviously, a follow-up sur- 
vey simultaneous to the blind survey could not be done after 
the new optimised survey reduction was completed. 

Using this method the follow-up sample was found to be 
78% complete above 50 mjy, 86% complete above 70 mjy, 
and 93% complete above 100 mjy over the whole survey 
area. Note that, as shown by the map in Figure [T] there are 
some regions with up to 95% completeness down to 50 mjy 
and 100% completeness above 70 mjy. 

3.2 Spectra 

As already noted in § 12.11 no 4.8 and 8.6 GHz follow-up 
observations have been made above 5 = —15°. The 3538 
sources at 5 < —15° followed up at 20 GHz were also ob- 
served at 4.8 and 8.6 GHz within a couple of weeks and 3332 
have good quality flux density measurements. The almost 
simultaneous observations ensure that the observed spec- 
tral behaviour is not significantly affected by variability. The 
AT20G pilot survey observations showed that, ona~ 1 yr 
timescale, the median variability index of sources selected 
at ~ 20 GHz is 6.9%, and only ~ 5% of sources vary their 
total flux density by more than 30% (Sadler et al. 2006). 

All but 27 AT20G sources have a counterpart in the 
NVSS (Condon et al. 1998), SUMSS (Mauch et al. 2007), or 
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Table 1 . Euclidean normalized differential source counts (log5 5 / 2 (iAr/(iS[Jy 1 /' 5 sr- 1 ]) for the AT20G sample as a whole and for the steep- 
(a < —0.5) and flat-spectrum (a > —0.5) sub-populations. Since spectral indices are different in different frequency ranges we need to 
specify the frequencies used for the source classification. Since spectral information is not complete, counts of each sub-population were 
corrected multiplying them by the factor iVtot /(-^flat + ■Wsteep) (see text). 
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MGPS-2 (Murphy et al. 2007) catalogues at ~ 1 GHz. The 
27 sources without a counterpart are discussed in Murphy 
et al. (2010) and in §|4] Combining our measurements with 
the ~ 1 GHz data we computed spectral indices, off (S oc 
v a \ between different pairs of frequencies (see §[4]), and used 
them to subdivide the sample into flat- and steep-spectrum 
sources adopting the usual boundary value a"^ = —0.5 (Wall 
1977). Although the measurements at ~ 1 GHz are far from 
simultaneous with those at higher frequencies, variability is 
not such a big problem at 1 GHz because the flux densities 
will be dominated by the stable extended steep-spectrum 
component of the source. 

The main effect of variability is the bias in favour of 
sources in a bright phase at the selection frequency of 20 
GHz. This will make some of the variable (and mostly flat- 
spectrum ) sources even more inverted but it would not move 
many sources from the steep-to the fiat-spectrum category. 

The error in the spectral index between pairs of fre- 
quencies depends on the flux density error and the fre- 
quency range. For the weakest sources typical rms errors 
are Sal = 0.14, Saf = 0.09, 5af = 0.06. For a\ the for- 
mal noise errors are smaller but the actual errors may be 
dominated by variability and are hard to quantify. In the 
following analysis we use all these ranges since the observed 
spectral curvature is significant. 

3.3 Source counts 

Figure [2] shows the Euclidean normalized total counts, and 
also the counts of flat- and steep-spectrum sub-populations. 
The latter counts are corrected for the incomplete spectral 
information. Since the incompleteness is due to observing 
constraints, unrelated to source properties, we correct it 
by simply multiplying the counts of each sub-population 
by JVtot /(iVkat + A^tcep). The spectral classification is based 



on the 1 to 5 GHz spectral index in the upper panel, and 
on the 8 to 20 GHz spectral index in the lower panel. The 
counts of steep-spectrum sources above ~ 3 Jy are unde- 
fined because of the lack of spectral information for extended 
sources (see § 12. 2[) . There is a clear increase of the steep- 
spectrum source counts based on a|°, compared to those 
based on af, as a result of spectral curvature steepening 
of high-frequency spectra. The lines show the predictions 
of the De Zotti et al. (2005) model for flat-spectrum (dot- 
ted) and steep-spectrum (dashed) sources, and for the total 
counts. The consistency with the observed total counts is 
remarkably good down to ~ 0.3 Jy. At fainter flux densi- 
ties the model shows a faster than observed convergence of 
the counts of flat-spectrum sources, and, as a consequence, 
under-predicts the total counts. The model counts of steep- 
spectrum sources are close to the observed ones (except at 
the brightest flux densities, where the measured flux den- 
sities are affected by resolution effects) if the spectral clas- 
sification relies on cvf, but are below those based on a|°. 
This is not surprising, given that the simple power-law spec- 
tra adopted in the model are not adequate to describe the 
complexity of real spectra. On the other hand, Vieira et al. 
(2010) find that their 150 GHz selected synchrotron emitting 
radio sources are consistent with a flat spectral behaviour 
(or q ~ —0.1) between 5 GHz and 150 GHz, and indeed 
the model is found to provide a very good fit to their 150 
and 220 GHz counts of radio sources. These apparently con- 
flicting results may be reconciled by noting that a spread 
of spectral indices causes the effective spectral index, a e ff, 
connecting the counts at the frequency v\ to those at the 
frequency U2, to shift with frequency. In the case of a Gaus- 
sian distribution of spectral indices with dispersion a a and 
of power-law counts with differential slope f3, the variation 
of the effective spectral index is given by (Kellermann 1964; 
Danese & De Zotti 1984; Condon 1984): 
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Figure 3. Euclidean normalized differential source counts for 
the 20 GHz AT20G sample (asterisks), for the WMAP 23 GHz 
NEWPS sample (squares, Massardi et al. 2009), for the 15.2 GHz 
9C sample (diamonds, Waldram et al. 2003, 2009), and for 8.4(x) 
and 4.8 (+) GHz surveys (the two latter have been multiplied 
respectively by a factor 10 and 100 to avoid overcrowding; see De 
Zotti et al. 2010 for references). Note that none of the counts have 
been re-scaled to match the AT20G survey frequency. The solid 
line shows the predictions of the De Zotti et al. (2005) model at 
20 GHz, while the dotted lines are, from bottom to top, the pre- 
dictions by De Zotti et al. (2005) at 15.2 and 8.4, and by Massardi 
et al. (2010) at 4.8 GHz. 



Figure 2. Euclidean normalized differential source counts for the 
whole AT20G sample (diamonds), for fiat- (asterisks) and steep- 
spectrum (triangles) sub-populations. The spectral classification 
is based on a\ or Og in the upper or lower panel, respectively. 
The counts of each sub-population were multiplied by the fac- 
tor iVtot / ( A'flat + iVsteep ) to correct for the incomplete spectral 
information (see text). The lines show the predictions of the De 
Zotti et al. (2005) model for the whole population (solid line), and 
for flat- and steep-spectrum sub-populations (dotted and dashed 
lines respectively) 



Aa e ff = 03 - 1) * (Ta *ln(l/ 2 /^l)- (1) 

For fiat-spectrum (a > —0.5) AT20G sources we find, for 
example, a a = 0.305 between 5 and 20 GHz so that Aa e ff is 
~ 0.19 for an Euclidean slope, /3 = 2.5. For comparison, the 
observed mean spectral indices steepen by (af ) — (af) = 
—0.15. This indicates that the effect of the spread of spec- 
tral indices may compensate for the steepening of the mean 
spectral index. Note that, strictly speaking, Aa e £f should be 
compared with observations for a sample selected at 5 GHz, 
not at 20 GHz. The latter selection is obviously biased in 
favour of 20 GHz bright sources, and therefore the value of 
(a|°) — (af) is bound to be larger (i.e. the observed steep- 
ening is bound to be smaller) than in the case of a 5 GHz 
selected sample. For the same values of a a and /3, Aa e ff be- 
tween 5 and 150 GHz is ~ 0.47. A more direct comparison 
between our results on source spectra and those by Vieira 
et al. (2010) is presented in § |4~T1 



The De Zotti et al. (2005) model is also in good agree- 
ment with the 31 GHz counts by Muchovej et al. (2009) at 
flux densities around 10 mjy, i.e. fainter than those reached 
by the AT20G survey, although at still fainter flux density 
levels it under-estimates the counts. As shown by Figure [2] 
the model predicts that, at flux density levels of a few tens 
of mjy, counts are dominated by steep-spectrum sources, 
consistent with the findings of Muchovej et al. (2009) . 

Figure compares Euclidean normalized differential 
source counts for the AT20G sample, for the WMAP 23 
GHz NEWPS sample (Massardi et al. 2009), for the 15.2 
GHz 9C sample (Waldram et al. 2003,2009) with lower (8.4 
GHz and 4.8 GHz) frequency surveys (see De Zotti et al. 
2010 for references). The predictions of the De Zotti et al. 
(2005) and Massardi et al. (2010) models have been overlaid 
for comparison. At higher frequencies the peak at high flux 
density due to the population of flat-spectrum AGNs is very 
evident and is confirmed by these observations. 

3.4 Compactness vs spectral properties 

Figure U shows the compactness estimated by the 6 km 
antenna visibility (see § 12. 3|) . versus a{ (upper panel) or 
q|° (lower panel). A Gaussian source of 0.3 arcsec in di- 
ameter has a visibility of 0.5 on a 4.5 km baseline at 20 
GHz. The unresolved sources on this plot (visibility > 0.85) 
are mostly much smaller than 0.1 arcsec so the range in 
distances in the sample has not smeared out this indica- 
tor of intrinsic size. 2793 of these sources also have good 
quality flux density measurements at the lower frequen- 
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Figure 4. Compactness parameter versus (upper panel) or 
o|° (lower panel). The compactness is quantified by the ratio 
of the 20 GHz flux density measured with the longest baselines 
(~ 4.5 km) to that obtained on the short baselines. Sources above 
a flux density ratio of 0.85 (i.e. above the horizontal dashed line) 
are unresolved even with the longest baseline and therefore have 
sizes smaller than about 0.2 arcsec. The vertical dashed line indi- 
cates a = —0.5, used as threshold to distinguish flat- and steep- 
spectrum sources. 



cies. There is a remarkably clean separation between steep- 
spectrum («i < —0.5) extended sources and flat-spectrum 
compact sources. Only 3% of the extended sources are flat- 
spectrum. Some of the them could be unrecognized Galac- 
tic sources with thermal free-free emission, but most are 
composite sources with the spectrum dominated by a flat- 
spectrum core but with enough emission in a jet or lobe 
to look extended. This class which includes upturning spec- 
trum sources will be discussed further in §[3] From Figure [3] 
the time-honoured choice of af = —0.5 as the boundary be- 
tween steep-spectrum/extended and flat-spectrum/compact 
radio sources, originally stemming from tentative indications 
of a double-peaked redshift distribution in the first survey at 
frequency > 1.4 GHz (Wall 1977), asserts its solid physical 
basis with the spectral index acting as a proxy for source 
size. 

The lower panel of Figure [4] shows the source com- 
pactness versus the 8-20 GHz spectral index. The compact 
sources show very similar behaviour but with a shift of the 
average spectral index to slightly steeper values and a tail 
extending to a steeper spectral index. The extended sources 
now show a much broader range in spectral index which is a 
combination of high frequency spectral steepening of the ex- 
tended sources and an increased number of hybrid sources 
with upturning spectrum due to extended sources with a 
flat-spectrum core which dominates at high frequency. 

It is also interesting that in both plots there is an appar- 
ent continuity of the compact sources across the -0.5 spec- 
tral index boundary. This strongly suggests that there is a 
population of steeper spectrum compact components (which 
includes the CSS sources) that are the tail of the compact 
flat-spectrum population. This is even more obvious in the 
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Figure 5. Radio colour-colour plot for the 3332 extragalactic 
sources with near simultaneous observations. Sources with single 
power-law spectrum lie on the dashed diagonal line. The central 
area of the diagram collects the flat-spectrum sources. The me- 
dian values have been overplotted for the extended sources (filled 
black diamond, see § l2.2l for details), the not-compact (filled black 
triangle) and compact (filled black square) sources classified ac- 
cording to the 6 km visibilities. The median values of spectral in- 
dices for bright (520GHz > 500 mjy, +) and faint (S20GHZ < 100 
mjy, X) are also plotted. 



lower panel, suggesting that a number of compact sources 
have spectral steepening at 20 GHz as is also seen in the 
95 GHz observations (Sadler et al 2008). 12% of the 20 GHz 
survey sample are compact and steep, according to this clas- 
sification, so while it's a significant extension of the popula- 
tion it has only a small impact on the models which divide 
the populations into flat- and steep-spectrum with the cut 
at -0.5. 



4 ANALYSIS OF TOTAL INTENSITY 
SPECTRA 

4.1 Distribution of spectral indices 

Figure [5] shows the color-color plot (Kesteven et al. 1977) 
comparing spectral indices at low and high frequencies. 3538 
extragalactic sources, the sub-sample with near simultane- 
ous data, were used in this analysis. 3332 sources with good 
flux density estimation were plotted in the Figure. 

It is clear from Figure [5] and from inspection of indi- 
vidual spectra, that a large fraction of the AT20G sources 
do not have power law spectra. We investigate the effect of 
this spectral curvature and the way it changes with flux den- 
sity in Table [5] which summarizes the percentage of sources 
with different spectral behaviour in various flux density se- 
lected sub-samples. The change of population composition 
with varying flux density is clearly visible: the flat-spectrum 
population is far more dominant at the highest flux densi- 
ties. 
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Table 2. Distribution of spectral behaviour for the 3332 sources with almost simultaneous 5, 8 and 20 GHz data and for the sub-samples 
with S20GHZ < 100 mjy, 100 ^ S20GHZ < 500 mjy, and S20GHZ > 500 mjy. See ig]for details. 



Spectrum 



Whole sample 
No. (%) 



< 100 mjy 
No. (%) 



100-500 mjy 
No. (%) 



> 500 mjy 
No. (%) 



Any 



3332 (100%) 1544 (46.4%) 1534 (46.0%) 254 (7.6%) 



> 0, al° > 



21) 

8 
v20 



< 



a 5 < °> a S > 
a% < 0, of < 



-0.5 < a§° < 0.5 



Inverted (I) 
Peaked (P) 
Upturning (U) 
Steep (S) 

Flat (F) 



195 (5.8%) 
183 (5.5%) 
102 (3.1%) 
1086 (32.6%) 



66 (4.3%) 
92 (5.9%) 
73 (4.7%) 
619 (40.1%) 



100 (6.5%) 
70 (4.6%) 
28 (1.8%) 
437 (28.5%) 



29 (11.4%) 
21 (8.3%) 
1 (0.4%) 

30 (11.8%) 



1766 (53.0%) 694 (45.0%) 899 (58.6%) 173 (68.1%) 



y 20 
„20 



-0.5 
-0.5 



Steep 
Flat 



1033 (31.0%) 
2299 (69.0%) 



578 (37.4%) 
966 (62.6%) 



422 (27.5%) 
1112 (72.5%) 



33 (13.0%) 
221 (87.0%) 



"5 < 
at > 



-0.5 
-0.5 



Steep 
Flat 



888 (26.7%) 
2444 (73.3%) 



543 (35.2%) 
1001 (64.8%) 



324 (21.1%) 
1210 (78.9%) 



21 (8.3%) 
233 (91.7%) 



a\ < -0.5 1 
af > -0.5 1 



Steep 
Flat 



876 (26.7%) 
2408 (73.8%) 



537 (35.3%) 
983 (64.7%) 



322 (21.3%) 
1191 (78.7%) 



17 (6.8%) 
234 (93.2%) 



1 Observations at 1 and 5 GHz are not simultaneous 
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Figure 6. Median spectral indices a|° in several flux density bins for some spectral populations, selected combining a| and a| with the 
same criteria as in Table [2] (see the text for details). The y-axis error bars have been calculated as errors on the median for a Gaussian 
distribution (i.e. 1.25/VN). 



It is also clear that the ratio of the flat- to steep- 
spectrum sources slightly decreases with increasing fre- 
quency, indicating a mild overall steepening of the spectral 
index of the population. The median spectral index over 
the whole sample is ag° = —0.28. The comparison with the 



median a| = —0.16 is a further indication of an overall 
steepening of the sample as the frequency increases. 

We can follow this behaviour to even higher frequen- 
cies for the synchrotron emitting sources detected in the 
150 GHz SPT survey. A cross-correlation of the AT20G sam- 
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pie with the catalog by Vieira et al. (2010) yielded 36 com- 
mon sources with S/N > 5 at both 20 and 150 GHz. Their 
median flux densities are 96 mjy at 20 GHz and 30.3 mjy 
at 150 GHz. Their median spectral index a|° = —0.12 is 
somewhat flatter than that for the whole sample, as an ef- 
fect of adding the 150 GHz selection. The mean spectral 
indices between 5 and 150 GHz and between 20 and 150 
GHz (al 5 = -0.38 and = -0.54) are substantially 
steeper than the values reported by Vieira et al. (2010): 
a\ 50 = —0.13 for the 57 sources for which they have ^ 5a 
detections at both frequencies, and CV20 = —0.31 for the 
top ten 150 GHz selected sources. However the mean spec- 
tral indices flatten to values close to those found by Vieira 
et al. (al 50 = -0.17 and a^,° = -0.30) if we empha- 
size the 150 GHz selection by taking the sub-sample with 
Si50GHz ^ 50mJy. We then conclude that the source spec- 
tral properties are strongly affected by the selection fre- 
quency. AT20G sources, selected at 20 GHz, mildly steepen 
between 5 and 20 GHz, while SPT sources selected at 150 
GHz exhibit a very flat spectral index across the frequency 
range 5-150 GHz. 

The only remaining significant difference with Vieira et 
al. concerns 0:5°; for the top ten 150 GHz selected sources 
they find a|° = 0.25, to be compared with a|° = —0.14 for 
the full AT20G/SPT sample, and to af = 0.03 for the sub- 
sample with 5150GHz ^ 50mJy. This difference is likely due 
to the use of different 5 GHz data: we used AT20G 5 GHz 
fluxes, while they used their own 5 GHz measurements, and 
the two sets of measurements were carried out at different 
epochs. 

By cross-matching the AT20G sample with NVSS, 
SUMSS, and MGPS-2, respectively at 1.4 and 0.834 GHz, 
we can extend the comparison to lower frequencies in Ta- 
ble H The median is a? = -0.16. However, the ~ 1 GHz 
observations are not coeval with the AT20G ones and for 
this reason spectral indices may be affected by variability, 
and the average spectral index may be biased towards high 
values because of the high-frequency selection. The compar- 
ison of the spectral indices over the three frequency ranges 
are plotted against the 20 GHz flux density in Figure [6] 
for the four spectral populations. This shows that the over- 
all spectral steepening is more prominent for the brightest 
sources (see also Figure [7] and [HJ that are typically flatter 
than the fainter sources (at any frequency): the fraction of 
flat-spectrum sources is 87 ± 6% if S20GHZ > 500 mjy, while 
it is 63 ± 2% if & GHz < 100 mjy. 

4.2 The different spectral types 

In Figure [6] we compare the change in spectral properties 
with flux density. The average difference in the median spec- 
tral index curves in this plot is a simple selection effect for 
the different spectral types, but the variation with flux den- 
sity is not and shows interesting trends. At high frequencies 
the spectrum of the flat-spectrum sources is essentially in- 
dependent of flux density while the spectrum of the steep- 
spectrum sources becomes steeper at lower flux densities 
at all the frequencies. This effect may be partly caused by 
some migration of sources between the flat- and the steep- 
spectrum source cut since the steep-spectrum sources are 
an increasing fraction of the low flux density sample. The 
sources with inverted and upturning spectra behave as ex- 
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Figure 7. Median values of of, and of (diamonds and black 
solid line, from top to bottom panel) for different flux density bins. 
The y-axis error bars have been calculated as in Figure [6] The 
dotted lines show the values of the first and third quartiles of the 
distributions in each flux density bin. 
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Figure 8. Median values of the spectral curvature o|° — a| for 
different flux density bins. 



pected if they are composite sources with a flat-spectrum 
core and steep-spectrum jets or lobes. At high frequency the 
core dominates and they have no dependence on flux density 
while their low frequency spectrum behaves like the steep- 
spectrum class of sources. Apart from the larger range of 
spectral differences from high to low frequencies the peaked 
sources have little changes of spectral indices with flux den- 
sity and look more like the flat-spectrum source population. 

The distribution of spectral indices is asymmetric with 
a tail towards steep spectra (see Figure [7} . The asymme- 
try is more clearly observable at the brightest flux densities 
and at the lower frequencies (where the difference of spectral 
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Table 3. UIS found in the AT20G sample with 5-8 GHz spectral index > +2.5. 



name 


■520GHz 


S&GHz 


•55GHz 


•5nvss 


a 20 
"s 


a 5 


err(a|) 


comments 




[mJy] 


[mJy] 


[mJy] 


[mJy] 










J070949-381152 


86+3 


33+3 


6±1 


< 2.5 


+1.1 


+3.0 


0.3 


UIS, 2MASS galaxy with K=12.7mag 


J111246-203932 


78±5 


30±2 


6±1 


< 2.0 


+1.1 


+2.9 


0.3 


UIS, No obvious optical ID 


J143608-153609 


113±6 


53±3 


10±1 


8.0 


+0.9 


+3.0 


0.2 


UIS, No obvious optical ID 
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Table 4. Matrix of the number of sources classified according 
to both the total intensity and polarisation spectral behaviour. 
The rows are the spectral shape in polarisation, the columns the 
spectral shape in total intensity. The spectral types are defined 
in Tabled 
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Figure 9. Percentage of flat- (asterisks) and steep- (triangle) 
spectrum sources defined between 8 and 20 GHz together with 
flat- (diamonds) and steep- (squares) spectrum sources defined 
between 1 and 5 GHz for different flux density bins. 



behaviour between faint and bright sources is more signifi- 
cant). This is a clear indication that extrapolations of flux 
densities from a low to a high frequency using only the mean 
or the median value for the spectral index are incorrect. A 
more correct approach should consider a flux density de- 
pendent spectral index distribution. Furthermore, extrapo- 
lations performed over large frequency intervals should also 
consider the spectral steepening which is significant at all 
flux density levels but most pronounced at high flux densi- 
ties (see Figure [HJ. 

The difference between the predicted and the measured 
ratio of flat- and steep-spectrum sources at the lower flux 
densities could be caused by the pair of frequencies used to 
make the separation. This would not change in the mod- 
els which assume only power law spectra. Figure [9] shows 
the fractions of flat- and steep-spectrum sources selected at 
the two extremes of our frequency range. There is a similar 
strong flux density dependence with an offset between the 
low and high frequency ranges. 



4.3 Ultra inverted spectrum sources 

Murphy et al. (2010) introduced the class of Ultra-Inverted 
Spectrum (UIS) radio sources which have a spectral index 
q§° > 0.7 and usually no detectable low frequency coun- 
terpart. In Fig [5] we see an even more extreme attribute of 
a subsample of these sources. 4 sources in the sample are 
listed with a 5-8 GHz spectral index a| > +2.5. Three of 
these with good quality flux density estimates are listed in 
Table [3] with their spectral index at 5-8 and 8-20 GHz and 



1 sigma errors. These sources are very weak at the lower 
frequencies so the errors are relatively large. 

All sources are just consistent with the maximum pos- 
sible synchrotron self absorption spectral index of +2.5 at 
the lower frequencies but are not consistent with the maxi- 
mum free-free absorption spectrum of +2.0. This steep self 
absorption spectrum suggests that these sources have been 
confined in a particularly dense environment, or are very 
young. 



5 POLARISATION PROPERTIES 

Polarisation was detected at 20 GHz for 768 sources, 467 of 
which also have simultaneous detections at 5 and/or 8 GHz. 
For each of these sources it is possible to compute the spec- 
trum of the polarised emission. We constructed the matrix in 
TableUby comparing the total intensity and the polarisation 
spectral behaviours for the sources with almost simultaneous 
total intensity observations and polarisation detection at 5, 
8 and 20 GHz. Although the spectral indices in total inten- 
sity and in polarisation are correlated (see Figure [10)) . for 
many sources the spectral shape in polarisation is substan- 
tially different from that in total intensity. Flat (F) sources, 
which are, on average, brighter, also have a higher polar- 
ized flux density (although not a higher polarisation degree) 
and this translates into a larger number of detections, but 
the spectral shape of the detected signal is not easily pre- 
dictable. The behaviour of the steeper spectrum power law 
sources is more predictable. Figure [TU] indicates that sources 
with power law spectra in total intensity tend to preserve 
the same spectral shape in polarisation, whereas flat spectra 
in total intensity may correspond to a large range of possible 
spectral shapes in polarisation. 

Low frequency depolarisation will move sources above 
the diagonal in the matrix but this seems to be a minor effect 
in the 5-20 GHz frequency range given the large number 
of sources that have steep-spectrum in both total intensity 
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Figure 10. Comparison of 8-20 GHz spectral indices in polar- 
isation (y-axis) and in total intensity (x-axis). Filled and open 
larger symbols are median values as described in Figure [5] 



and in polarisation. Low frequency depolarisation is seen as 
a small excess of sources above the power law line in the 
steep-spectrum region of Figure [TO] 

The half of the matrix below the diagonal in Figure 
1101 includes a class of sources that show decreased polarised 
flux density at the higher frequencies and includes a distinct 
group of sources with positive spectral index in total in- 
tensity and negative spectral index in polarisation. These 
sources may be composites with a steeper spectrum po- 
larised component and a less polarised flat-spectrum core 
that dominates at higher frequency. 

To calculate the median fractional polarisation at 20 
GHz we have exploited the Survival Analysis techniques 
and the Kaplan-Meyer estimator (as implemented in the 
ASURV code, LaValley et al. 1992). The large number of 
upper limits (83 percent of the sample) at 20 GHz makes 
any Survival Analysis over the whole sample unreliable. To 
get meaningful results we need to restrict ourselves to the 
1059 sources with S20GHZ > 250 mjy, for which the number 
of upper limits (526) is < 50 percent of the sample. If only 
the robust detections are used the median fractional polar- 
isation degree is 2.6 percent. For these sources the Survival 
Analysis yields a mean polarisation degree of 2.7 percent, in 
perfect agreement with the mean polarisation degree at 22 
GHz measured by Jackson et al. (2010) for a complete sam- 
ple of extragalactic sources stronger than 1 Jy in the 5-year 
WMAP catalogue (2.7 ± 0.2 percent) and with the median 
polarisation degree at 18.5 GHz of southern flat-spectrum 
sources in the Kuhr et al. (1981) 1 Jy sample, found to be, 
again, 2.7 percent (Ricci et al. 2004b). 

There is no evidence of a higher fraction of integrated 
polarised flux density at low flux density levels. 

We can also use the compactness classification based on 
the long baselines to look at the polarisation characteristics 
of the compact and extended sources. This confirms the ob- 
servation that most steep-spectrum extended sources have 
the same spectral behaviour in both total and polarised in- 
tensity see Fig. I10|) . If we consider only the two groups of 




PA 8GHz [%] 

Figure 11. Comparison of polarisation angles at 5 and 8 GHz 
(upper panel) and at 8 and 20 GHz (lower panel). 



compact flat-spectrum and compact steep-spectrum sources 
the mean of the detected values is 2.9% and 3.8%, respec- 
tively. 

The ratio of fractional polarisation at different frequen- 
cies is a measure of depolarisation. The median of the ratio 
m^/m2o is 0.83 indicating on average a 17% depolarisation 
at 5 GHz at any flux density level. 

The polarisation angle is quite similar at all the frequen- 
cies (see Figure [TTjl . However, the position angle differences 
have errors which are too large to determine meaningful ro- 
tation measures and the scatter is consistent with the ex- 
pected intrinsic dispersion of rotation measure observed at 
lower frequencies. 

15 sources show polarisation fraction larger than 5 
percent with 20 GHz flux density larger than 1 Jy. The 
source with the brightest polarised flux density is the well- 
known blazar 3C279, with P20 =1.1 Jy and S20 GHz =20 
Jy. The source with the largest fractional polarisation is 
AT20GJ195817-550923 which is 19 per cent polarised with 
S20 GHz =0.058 Jy 
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6 CROSS-MATCHES IN OTHER SPECTRAL 
BANDS 

6.1 Lower radio frequencies 

Cross matches with catalogues at 1.4 GHz(NVSS) and 0.843 
GHz (SUMSS-MGPS-2) have been discussed in §g] and in 
Murphy et al. (2010). Such comparisons demonstrated that 
the dominant population changes with frequency and that 
extrapolations assuming power laws behaviour from GHz 
frequencies cannot reconstruct the 20 GHz population. 

Only 0.4% of the AT20G sources have no low frequency 
counterpart. They have non simultaneous spectral index be- 
tween 1 and 20 GHz larger than 0.7 and mostly belong to 
the inverted spectra sources group in the upper right quad- 
rant of Figure These belong to the class of sources with 
a|° > 0.7 that in Murphy et al. (2010) have been defined as 



'Ultra Inverted Spectra sources' 
belong to this class. 



1.2% of the AT20G sources 



6.2 Optical identifications and redshift 

Optical counterparts for the 4932 AT20G sources with |b| > 
10° were found by cross-matching the radio positions with 
optical positions in the SuperCOSMOS database (Hambly et 
al. 2001). Optical identifications were chosen to be the clos- 
est optical source to the radio position within a 2.5 arcsec 
radius and with a B magnitude brighter than 22, which is the 
SuperCOSMOS completeness limit (see Sadler et al. 2006). 
This resulted in 2958 (60%) sources with optical identifica- 
tions. The SuperCOSMOS classifications revealed that 1762 
(60%) are point-like and 1172 (40%) are extended (for the 
remaining sources no classification is available). The sim- 
ple position-based optical identification procedure is good 
for point sources but 2-3% of the sources with more com- 
plex radio structure are likely misidentified. A more detailed 
analysis of the optical identifications including source struc- 
ture will be made in a future paper (Mahony et al., in prepa- 
ration). 

The fraction of sources with optical identifications in 
the full AT20G sample is less than is seen in the AT20G 
Bright Source Sample (Massardi et al. 2008). This is due to 
the decrease in flux density limit from 500 mjy in the BSS to 
40 mjy for the full sample and corresponds to the increase in 
the number of galaxies as the 20 GHz flux density decreases. 

Redshift information for 825 (28% of those with optical 
identifications) sources was obtained by searching the NASA 
Extragalactic Database (NEE0). Of these, 443 are classified 
as point-like in the optical and 379 as extended. The median 
redshifts are 1.25 and 0.25 for each of these classes respec- 
tively. Similar fractions, 76% and 71%, of the point-like and 
extended sources respectively have flat spectra between 8 
and 20 GHz. 

Spectral index versus redshift is plotted in Figure [T2l for 
various flux density ranges. The most obvious effect is that 
most of the bright sources are flat-spectrum high redshift 
sources, while the closest sources are fainter and include a 
larger population of steep-spectrum sources. With redshifts 
for only 14% of the AT20G sources there are a number of 
strong selection effects. The strong flat-spectrum AGNs are 
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Figure 12. Spectral indices o|° versus redshift for different flux 
density intervals: S20GHZ > 0-5 Jy, diamonds; 0.5 > S20GHZ > 0-1 
Jy, triangles; S20GHZ < 0.1 Jy, circles. Filled symbols correspond 
to median values for each flux density selected sub-sample. 



bright quasars with easily measured redshifts. The galaxies 
in the sample are fainter and redshifts are strongly biased to 
the closer galaxies. For nearby sources (z < 0.2) the redshift 
surveys are providing redshifts and some complete subsam- 
ples are available. Sadler et al. (in preparation) are analysing 
the 203 AT20G sources in the 6dF Galaxy Survey. In the fol- 
lowing discussion we restrict our comments to redshift effects 
which are less influenced by sample completeness. 

Figure [T^] indicates that the flat-spectrum population 
shows a steepening of the a|° spectral index at high z but 
the steep-spectrum population does not show any signifi- 
cant change with redshift. As illustrated by Fig. 1131 how- 
ever, both populations show a slight increase with redshift 
of the spectral curvature at lower frequencies (q|° — a|). 
This behaviour is similar for all flux density ranges and is 
the expected effect of the K correction on a population of 
sources with curved spectra. The redshift distribution is sim- 
ilar for flat- and steep-spectrum sources. We have redshift 
measurements for 298 flat-spectrum and 147 steep-spectrum 
sources; the median values are 0.83 and 0.93, respectively. 

6.3 Cross-matches in the X-ray band 

The ROSAT All-Sky Survey Bright Source Catalogue 
(RASS-BSC, Voges et al. 1999) includes 8896 X-ray sources 
in the Southern sky with |fe| > 5°. 280 sources from the 
whole AT20G sample (4.6%) have a counterpart within 10 
arcsec. The median 20 GHz flux density for the X-ray de- 
tected sample is 189 mjy. The median a|° for the 161 ex- 
tragalactic sources with data on the three AT20G follow-up 
frequencies is -0.20; the median a% of the same sources is 
—0.06. Only 76 sources with flux density larger than 189 
mjy have a polarisation detection. The median fractional 
polarisation is 2.7%, equal to the average calculated over 
the whole sample. There is no clear relation between X-ray 
and radio flux densities in Figure [14] However, the X-ray 
detection rate for AT20G survey sources increases with in- 
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Figure 13. Variation of spectral index o|° — a| vs redshift. Sym- 
bols are the same as in the previous image. 
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Figure 15. Percentage of X-ray detections at varying 20 GHz flux 
densities for the whole sample (diamonds), and for flat- (asterisks) 
and steep-spectrum (triangles) sources. 



sources leaves 540 extragalactic 7-ray sources in the south- 
ern sky. Of these, 233 (43%) have an AT20G counterpart 
within the 95% confidence ellipse. This is 4% of the whole 
AT20G sample - very similar to the X-ray detection rate 
but only 56 (1%) of the AT20G sources have both X-ray 
and 7-ray detections. The detection rate also increases with 
7-ray flux density such that the brightest 7-ray sources are 
all detected in the AT20G survey. The median 20 GHz flux 
density for the 7-ray detected sources is 519 mjy, signifi- 
cantly brighter than the X-ray detections. Figure [TBI shows 
the spectral index distribution of the Fermi-AT20G detec- 
tions compared to the full AT20G sample. As expected, 
the sources that have gamma-ray counterparts are predom- 
inately fiat-spectrum sources with median spectral indices 
a|° = —0.16 and a| = —0.01. Kolmogorov-Smirnov tests 
confirm that these distributions are not drawn from the same 
population at the 99% confidence level. 



Figure 14. Comparison of 20 GHz flux densities with X-ray flux 
densities from the ROSAT-BSC (Voges et al. 1999). 



creasing 20 GHz flux density, as shown in Figure 1151 Here 
we binned the AT20G catalogue in flux density and plotted 
the percentage of those that were detected in the ROSAT 
dataset for each bin. There is a clear trend in the percent- 
age of detections, implying that the brightest sources have 
higher probabilities of showing a counterpart in the X-ray 
band and, as a consequence, a relation between the mecha- 
nisms that generate the flux densities in the two bands. 

6.4 Cross-matches with the 7-ray band 

These results are presented in more detail in Mahony et 
al. (2010) but we summarise the statistics here and make 
comparison with the X-ray detections. The first Fermi-LAT 
catalogue (1FGL) (Abdo et al. 2010) includes 1451 7-ray 
sources observed during the first 11 months of operation. 
Reducing this catalogue to the same region of sky covered 
by the AT20G survey and excluding known Galactic 7-ray 



7 CONCLUSIONS 

We analyzed a sample of 5808 extragalactic sources in the 
AT20G survey. 

We present the source counts for the whole sample and 
for spectrally selected sub-samples. As expected, our analy- 
sis has confirmed that the high frequency bright source sam- 
ple is dominated by flat-spectrum sources, while the fraction 
of steep-spectrum sources increases with decreasing flux den- 
sity. As a consequence, the mean spectral index steepens and 
the spectral index distribution becomes more symmetric at 
lower flux density. Extrapolations of the counts to different 
frequencies must consider flux density dependent steepening 
and changing shape of the spectral index distributions. 

There is an overall spectral curvature causing spectral 
steepening of the population with frequency, which is more 
prominent at higher flux densities and at higher frequencies. 
The steepening is stronger for the more compact sources. 
The spectral distribution of the population of flat-spectrum 
sources is not a power law and changes with both flux den- 
sity and redshift. Modelling this population will be a chal- 
lenge. 
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Figure 16. 5-20 GHz spectral index distribution for the full 
AT20G sample (top panel), those sources with Fermi counter- 
parts (middle panel) and those with X-ray detections (bottom 
panel). The dotted line at a = —0.5 denotes the commonly used 
dividing line to separate steep- and fiat-spectrum sources. 
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The median fractional polarisation at 20 GHz is 2.7% 
for 520 GHz > 250 mJy. Above this flux density limit the 
sample is dominated by flat-spectrum sources. We find no 
evidence of a variation of the polarisation degree with flux 
density but, because of the very large number of upper lim- 
its, we could not extend this analysis to lower flux densi- 
ties. Comparing the polarisation distributions for simulta- 
neous detections at 5 and 20 GHz we find an indication 
of an average 17% depolarisation at 5 GHz independent of 
flux density. Although spectral indices in total intensity and 
in polarisation are correlated, many, mainly flat-spectrum, 
sources have a spectral shape in polarisation substantially 
different from that in total intensity. However, sources with 
steep power law spectra in total intensity tend to preserve 
the same spectral shape in polarisation. 

189 sources have a counterpart in the ROSAT-BSC cat- 
alogue and there is a general trend of increasing X-ray de- 
tection rate of AT20G sources with increasing 20 GHz flux 
density. There is no strong correlation between X-ray and 
radio flux densities. The distribution of radio spectral in- 
dices for the X-ray detected sources is similar to that for 
the total AT20G sample. 

There are 233 AT20G sources that have gamma-ray 
counterparts in the first Fermi-LAT catalogue. The major- 
ity of these sources have flat spectral indices between 5 and 
20 GHz which is consistent with the high energy emission 
originating from compact cores of AGN. 

Although both X-ray and 7-ray detection rates are re- 
lated to the 20 GHz flux density, the 7-ray emission is more 
closely correlated with both the radio flux density and its 
spectrum. This suggests that the 7-ray and radio emission 
processes are more closely linked in space and time while the 
X-ray - radio link may be more indirect with AGN activity 
occuring in the same sources over a period of time. 
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